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ABSTRACT 

Context. Grain growth in circumstellar disks is expected to be the first step towards the formation of planetary systems. There is now evidence 
for grain growth in several disks around young stars. 

Aims. Radially resolved images of grain growth in circumstellar disks are believed to be a powerful tool to constrain the dust evolution models 
and the initial stage for the formation of planets. In this paper we attempt to provide these constraints for the disk surrounding the young star 
CQ Tau. This system was already suggested from previous studies to host a population of grains grown to large sizes. 

Methods. We present new high angular resolution (073 - 079) observations at wavelengths from 850 yum to 3.6 cm obtained at the SMA, 
IRAM-PdBI and NRAO-VLA interferometers. We perform a combined analysis of the spectral energy distribution and of the high-resolution 
images at different wavelengths using a model to describe the dust thermal emission from the circumstellar disk. We include a prescription for 
the gas emission from the inner regions of the system. 

Results. We detect the presence of evolved dust by constraining the disk averaged dust opacity coefficient ft (computed between 1.3 and 7 mm) 
to be 0.6±0.1. This confirms the earlier suggestions that the disk contains dust grains grown to significant sizes and puts this on firmer grounds 
by tightly constraining the gas contamination to the observed fluxes at mm-cm wavelengths. We report some evidence of radial variations in 
dust properties, but current resolution and sensitivity are still too low for definitive results. 
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1. Introduction 

In recent years the study of planetary system formation has 
grown enormously. While it is quite firm that the formation 
of planets should occur during the pre-main sequence (PMS) 
phase of the central star, when the circumstellar disk is still 
rich in gas and dust, our understanding of the complex variety 
of factors which contribute to that process is still very sketchy. 
Two possible scenarios are currently considered for planet for- 
mation: the "disk instability" for gas giant planets (see e.g. 
Durisen et al.|2007 1, and the "core accretion " for solid plane- 
tary cores from the dust particles in the disk (see e.g. Lissauer 



1993] >. The growth by coagulation of dust particles in the cir- 
cumstellar matter from the typical sub-micron size found in the 
interstellar matter (ISM) to millimeter and centimeter sizes is 
now considered as the first step in the latter scenario; constrain- 
ing observationally the dust evolution in the disk midplane al- 
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low us to put constraints on the planet formation theories (see 
|Beckwithet al.|2000) . 

Following the initial attempts to characterize the dust pop- 
ulation using the submm spectral energy distribution of disks 
(e.g. [Beckwith & Sargent|[T99l] |Miyake & Nakagawa|[l993l >, 
the last decade has seen a number of studies dedicated to the 
search for large grains in protoplanetary disks at mm and cm 
wavelengths. These studies have been made possible by the 
improved capabilities of millimeter and centimeter wave ar- 
rays, but are still limited to bright objects (Wilne r~et al.|2000 

hnArL Fr 7° T — l l ^rvr> 1 I l ^r>r>--> I FH — "i 7 — ^^T^^^R^^^ ^ l 



[2005l|Testret al.|200H[2003i|Rodmann et al.|2006HNatta et al.| 



|2007t |Lommen et al.||2009| |Ricci et aX]|20lO| l. Even in these 
favourable cases it has been so far impossible to resolve the 
radial variations of the dust properties and high angular resolu- 
tion observations have mostly been used to constrain the mass 
distribution in disks by resolving the surface brightness profile 
at a single wavelength (|Isella et al.|2007| [Hughes et al. 2008 ; 



And rews et al.||2009| |Isella e t al. 2009f7 Resolving the radial 



variations of the dust properties is particularly interesting to 
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constrain models of dust evolution in protoplanetary disks that 
include fragmentation and radial drift of large grains (see e.g. 
Brau er et al.|2008|[BFnstiel et al.|2010a) . Very recently peHa 
et al.| ( |2010| ) used high sensitivity and high angular resolution 
CARMA observations at 1.3 and 3 mm to set limits on the pos- 
sible radial variations of the dust properties in the protoplane- 
tary disks around RY Tau and DG Tau. 

In this paper we focus on the protoplanetary disks of CQ 
Tau in the Taurus-Auriga star-forming region. CQ Tau is one of 
the closest (100 pc, Hipparchos, Perryman et al.|1997| l young 
stellar objects with disk observable from the northern hemi- 
sphere, where most of the current high sensitivity and high 
spatial resolution millimeter and centimeter wave arrays are lo- 
cated. The central star is a well known variable of the UX Ori 
class with spectral type A8 and estimated age ~ 10 Myr. The 
star-disk system has been studied in detail by various authors 
and the presence of grain growth to large sizes (~cm) was in- 
ferred from the mm range of the SE P (|Chiang et aT]|2001| 
Natta et al.|2000}|Tesfi et al.|2001||2003"t|Chapillon et al.|2008) ' 
Here we refine the investigation of dust properties in the disk 
of CQ Tau using new high-resolution observations obtained 
from the VLA, the PdB and the SMA interferometers at cm 
and mm wavelengths (Sections [2] and [3]l. We combine our new 
observations with archival and literature ones and re-analyze 
them together. We use a dusty disk model with an accurate 
parametrization for the dust opacity. We include a prescription 
to account for the gas emission. The average grain properties 



are first constrained by fitting the observed SED (Section 4.3 I, 
separate model fits to the visibilities at various wavelengths are 
then used to contrain the disk structure and to probe possible 
radial variations of dust properties (Section[5]l. 

2. Observations and data reduction 



2.1. VLA observations 

The new observations at 1.3 cm (22.4 GHz, K-band) and 3.6 cm 
(8.4 GHz, X-band) were done on September 1 1 and 14, 2006 
and November 03, 08 and 12, 2007 with the NRACQVLA in 
Soccorro. The array was in the B configuration, with a range of 
baselines from 210m to 11.4 km providing a maximum angular 
resolution of 0.7 and 0.3 arcsec at 3.6 and 1.3 cm respectively. 
CQ Tau had been already observed in both bands with the VLA 
at lower angular resolution, but only an upper limit at 3.6 cm 
was available from previous observations (Testi et al. 2001| l, 
while the previous 1.3 cm observations were not published. To 
maximize the signal-to-noise we decided to combine all VLA 
data (see Table[T|. The data reduction was performed using the 
CASA software package, following the standard procedures as 
described in the Data Analysis Cookbool^] At the time of this 
work CASA is still under development so we compared step 
by step the reduction products using AIPS and with our work 
we confirm that the two packages yield to consistent results. 



1 The National Radio Astronomy Observatory is a facility of the 
National Science Foundation operated under cooperative agreement 
by Associated Universities, INC. 

2 The Common Astronomy Software Applications (CASA) manual 
is available on casa.nrao.edu 



Standard calibration procedures were employed, the complex 
gain calibrators used and their fluxes as derived from compari- 
son with the VLA primary calibrators observations are reported 
in Table[T] flux calibration is expected to be accurate within 5% 
at 3.6 cm and 10% at 1.3 cm. For deconvolution and imaging 
the natural weighting of visibilities was used in all cases. Once 
calibrated, all datasets belonging to the same band were added 
together for the final deconvolution and imaging, so that the 
final maps of CQ Tau at 1.3 and 3.6 cm result from the combi- 
nation of all VLA data calibrated in CASA (see PigJTJ. 

2.2. PdBI observations 

Observations of CQ Tau were carried out simultaneously at 
2.7 mm and 1.3 mm in the most extended configuration of the 
six-element IRAM Plateau de Bure interferometer on Jan 15, 
2006. Conditions during the observations were excellent as tes- 
tified by the high atmospheric transparency (pwv= 0.5-1 .0 mm 
for 70% of the time) and high atmospheric phase stability 
(o"0 < 40° at 1.3 mm) observed over the 7.5 hours track. The 
spectral configuration was set to cover a contiguous band of 
550 MHz for maximum continuum sensitivity. The RF calibra- 
tion was measured on 3C454.3, and amplitude and phase cali- 
bration were made on 0528+134. The absolute flux calibration 
was determined using as primary calibrator MWC349. Data re- 
duction and calibration were made using the GILDAS software 
package in the standard antenna based mode. Continuum point 
source sensitivities of 0.15 mJy/beam and 0.5 mJy/beam were 
obtained at 2.7 mm and 1.3 mm, respectively, and consistent 
with the measured system temperatures (250-300 K). 

2.3. SMA observations 

The CQ Tau system was observed with the Submillimeter 
Array (SMA) on Mauna Kea, Hawaii, on 2008 Jan 24. The 
array was in the extended configuration, and the six available 
antennas provided baseline lengths from 60 to 225 meters. The 
atmospheric transparency was excellent, with 225 GHz opac- 
ity of 0.05 as measured at the nearby Caltech Submillimeter 
Observatory. The double sideband receivers were tuned with a 
local oscillator (LO) frequency of 340.8 GHz, and the correla- 
tor provided 2 GHz of bandwidth in sidebands +5 GHz from 
the LO frequency. Useable data were obtained for CQ Tau over 
an hour angle range of -0.9 to +2.2. Complex gain calibra- 
tion used observations of the calibrators J0530+135 (1.75 Jy) 
and 3C1 1 1 (3.0 Jy) interleaved with CQ Tau. The passband re- 
sponse was calibrated using observations of the strong source 
3C273. The flux scale was determined by bootstrapping obser- 
vations of 3C273 calibrated against planets and is accurate at 
the 15% level. All of the data editing and calibration were per- 
formed using standard routines in the MIR software package 
for IDL. 

3. Results 

In Figure [T] we show the interferometric images of the CQ Tau 
system from 0.87 mm through 3.6 cm. We detect strong emis- 
sion at millimeter wavelengths and the emission is resolved at 
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Table 2. Summary of available data at mm-cm wavelengths; 
the new observations presented in this work are highlighted. 
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a (J2000.0) Epoch 2000.00 
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Fig. 1. CQ Tau continuum maps in mm-cm bands from our new 
high angular resolution data (from top left to bottom right): 
0.87 mm, 1.3 mm, 2.7 mm, 7 mm (from lTesti et al.|2 003 ), 1.3 
cm, 3.6 cm maps. In each plot the position of the central star 
and the beam ellipse are shown. Contours levels are spaced by 
2cr where <x is 9 mJy at 0.87 mm, 1 mJy at 1.3 mm, 0.3 mJy at 
2.7 mm, 0.1 mJy at 7 mm, 0.02 mJy at 1.3 cm and 0.009 mJy at 
3.6 cm. All maps have been corrected from their original epoch 
to the epoch 2000.0, using the proper motion of the optical star 
as measured by Hipparcos (RA: 2.67, Dec: -26.17 mas/year), 
the star is also plotted at the epoch 2000.0 coordinates, as given 
by Simbadat CDS. 



all wavelengths from 0.87 mm through 7 mm (as discussed in 
|Testi et al.|2003l >. 

At 1.3 mm, we measure a spatially integrated flux of 103 
mJy, which is ~40% smaller than that measured at 2"resolution 
with the same instrument (IRAM-PdBI, |Chapillon et al.|2008) . 
Even if this discrepancy may be marginally within the cali- 
bration uncertainties at this wavelength (10-20%), we caution 
that the difference may be caused by the presence of extended 
(> 3") dust emission partially filtered by our observations. This 
however does not seem to be consistent with measurements at 
the other wavelengths (which do not show extended emission) 
and our model fit to the SED (see Sect. |4.3| l which does not 
seem to require additional extended components. 



A 

(mm) 


Synthesized beam|^] 
(arcsec, deg ) 


Pealj*] 

(mJy) 


Flux 

(mJy) 


rrnij^] 

(mJy) 


Ref[] 


0.87 


0.85 x 0.45, 95 


253 


421 


9 


(1) 


1.3 






221 


40 


(2) 


1.3 


1.5 




162 


2 


(3) 


1.3 


2 




143 


8.4 


(4) 


1.3 


0.45 x 0.29, 35 


18 


103 


1 


(1) 


2.7 


0.91 x 0.56, 34 


9 


22 


0.3 


(1) 


2.7 






18 


0.6 


(2) 


3.4 






13.1 


0.5 


(3) 


6.9 


0.92 x 0.72, -37 


1.0 


2.0 


0.1 


(5) 


13.4 


0.33 x 0.30, 2 


0.08 


0.28 


0.02 


(1) 


35.7 


0.91 x 0.81, -15 


0.046 


0.079 


0.009 


(1) 


62.5 


31 x 23,49 


<0.75 




0.25 


(6) 



" The first two values correspond to the size of the axes of the syn- 
thesized beam ellipse, while the third to its position angle. The beam 
is reported only for interferometric observations 

* Only for interferometric observations 

c In the case of data reduced in this work the rms is extracted di- 
rectly from the residual noise in final cleaned image, otherwise it is 



d Data from the literature: (1): this work, (2): Mannings & Sargent 


jl997i, (3): Chapillon et al. 


(20081, (4): Mannings & Sargent (2000), 


(5): Testi et al. (2003 1, (6):" 


festi etal. 1 2001 1. 



At centimeter wavelengths, our new VLA observations al- 
lowed us to detect the system for the first time at wavelengths 
longer than 7 mm, with a detection at both 1.3 and 3.6 cm. The 
source is not resolved at 3.6 cm; at 1.3 cm we tentatively de- 
tect extended emission at a ~2cr level around the optical star 
position. In Table [2] we report the beam sizes and sensitivities 
of our maps as well as peak and integrated fluxes for CQ Tau, 
in the same table we also report the fluxes from the literature 
that we have used in Sect. l4.3l for the fits to the SED. 

The possible detection of an extended ring-like emission 
at 1.3 cm coincident with the disk detected at shorter wave- 
lengths, if confirmed, suggests that at this wavelength we are 
still detecting the dust thermal emission. The 3.6 cm emission 
is unresolved and it is consistent with emission from ionized 
gas in the inner region of a wind or around the central star. 

After adjusting the astrometric calibration of the radio 
observations to account for the measured proper motion of 
CQ Tau from Hipparchos (Perryman et al. 1997), the optical 
position of the star is coincident with the peak emission at (sub- 
millimeter wavelengths. 

4. Dusty disk model, gas emission and SED 

4.1. Dusty disk 

To analyze the data we have used the model for passive flared 
two-layer disks presented by Dullemond et al. ( 2001 j >, based on 
|Chiang & Goldreich ( 1997 1. The model requires as input a set 
of parameters for the central star and the disk, which are used 
to compute the disk structure, SED and radial intensity profile 
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Table 1. Summary of PdBI and VLA data (new + archival) reduced for this work. For each dataset we show the observation 
date, the array configuration, the primary flux calibrator used and the flux assumed from the observatory calibration plan, the 
amplitude and phase calibrator used and its derived flux at the time of the observations (note that we reported the calibrator 
names as in the databases of the observatories at the time of the observations, 0528+134 used at PdBI in 2006 is the same 
object as 0530+135 used at the VLA in 1999). 



Band 


Project 


Date 


array 


Flux calibrator 


Gain calibrator 








conf. 


Name 


Flux (Jy) 


Name 


Flux (Jy) 


1.3 mm 


P054 


06 Jan 15 


A 


MWC349 


1.9 


0528+134 


2.1 


2.7 mm 


P054 


06 Jan 15 


A 


MWC349 


1.2 


0528+134 


3.1 


1.3 cm 


AK0475 


99 Jan 14 


C 


0137+331 


1.117 


0530+135 


2.77 




AT0327 


06 Sep 14 


B 


0137+331 


1.117 


05595+23539 


0.292 




AT0352 


07 Nov 03 


B 


0137+331 


1.117 


05595+23538 


0.246 


3.6 cm 


AT0327 


06 Sep 11 


B 


0137+331 


3.164 


05595+23539 


0.463 




AT0352 


07 Nov 08 


B 


0137+331 


3.164 


0547+273 


0.242 




AT0352 


07 Nov 12 


B 






0547+273 


0.24^] 



" In this dataset the flux calibrator is absent; for calibration we assumed the flux value of the phase-calibrator derived in the same band 
four days before (07 Nov 08). 




Fig. 2. Results of the SED fitting using the reference dust described in the text (Section|4j). Left: confidence limits on dust opacity 
index j3 and gas emission index ag; contours represent the 68% and 95% confidence levels for the fit (i.e. they represent the 
Ax 2 = 2.3, and 6.17 values). Middle and right: examples of fits to the SED within the 68% confidence level limits; the dusty disk 
model is displayed as a continuum line, the free-free power law for the gas emission as a dotted-dashed line, and the total model 
(dusty disk + free-free) as a dashed line, the flux upper limit at 6 cm as an arrow. The middle and right panels are for q-jj = -0.2 
and 2.0, respectively. 



at each wavelength. The stellar properties are assumed known 
as in |Testietal.| ( |2003] l (T eff =6900 K, L*=6.6 L , M* = 1.5 Mq, 
D=100 pc). The disk surface density distribution is assumed as 
£(R) = 2/ioau (R/40AU)~'\ where X40AU is the surface density 
at the fiducial radius of 40 AU. For the dust opacity at a given 
frequency k y we consider the parameterization by Miyake & 
as: k v — ( J «(a)fl 3 k Vt i(a)da) / (J n(a)a i da), 



Nakagawa 



(1993 



where a is the particle size distributed between a minimum and 
a maximum value a min and a max according to n(a) oc a~ q , and 
ky \ is the opacity per unit mass of the single particle, which 
depends on grain shape and chemical composition. The value 
of a,,,,,, is set to 0.01 fim according to the typical value observed 
for unprocessed dust in the ISM, while for the power law ex- 
ponent we have used q = 3 as proposed in previous studies 



and theoretical expectations (see|Draine 2006 1. In this work we 



assume as reference the same ensemble of dust used in IRiccil 
( 2010) 1, where porous spherical grains are composed by 



et al. 



(volume percentage) 7% silicates, 21% carbonaceous materi- 
als, 42% water ice and 30% vacuum ( |Pollack et al.|19 94). It is 
well known that the dust opacity, and thus the exact value of 
a max derived from fitting the millimeter dust emission, depends 
strongly on the dust composition and the power law index of the 
dust size distribution (e.g. |Natta et al.|2007||Ricci et al.|2010| >. 
In Sect. 4.3 we will show which is the range of uncertainty 
in our derivation of a max by exploring its dependence on the 
assumed dust model and the grain chemical/physical proper- 
ties: removing the water ice, varying the vacuum volume from 
10% to 50% and the grain size distribution index q from 2.5 
to 3.5. We stress that the observations constrain very well the 
dust opacity index in the simplified first-order approximation 
k v oc yP. In this work we compute the slope of the dust opacity 
/J from the ratio of the values of k v at 1.3 and 7 mm. This pro- 
vides a parameter that can be directly compared to other stud- 
ies, but note that the shape of k v that we use in the disk model 
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is the parameterization by Miyake & Nakagawa ( 1993) written 
above. The dependence of /3 on a max and q was discussed in 
|Natta & Testi] ( p004] >, and |Ricci etaE] pOTO] ) showed the case 
of our reference ensemble of dust. 



4.2. Gas emission 

Young stellar objects are known to show radio continuum emis- 
sion from circumstellar gas, this emission can either be ther- 
mal or non-thermal (see e.g. Andre et"aL|l987| Skinner et al. 
|1993| l. The radio emission can normally be well approximated 
by a power law of the form F A oc X" B with negative values 
of the radio spectral index (aff) normally attributed to thermal 
emission (Mezger & Henderson 1967; Panagia & Felli 1975; 



Franco et al. 2000 1 and positive values to non-thermal emis- 



sion ( |Guedel et al.|[l989] l. |Skinner et aT] ( fT993] l found spec- 
tral indices spanning the range between -1.2 and 1.2 in a sam- 
ple of Herbig Ae/Be stars. Using an equivalent parametrization 
Rodmann et al. ( 2006| l applied a simplified correction for free- 



free contamination at radio wavelengths (assuming a priori a 
50% free-free contribution at 1.3 cm), and found values for a$ 
between -0. 1 and 0.4 in a sample of low-mass PMS stars in the 
Taurus-Auriga star forming region. These results show the va- 
riety of spectral indices displayed by radio emission in young 
stellar objects and highlight the need of a proper characteriza- 
tion of the gas emission to accurately constrain the dust prop- 
erties from millimeter and centimeter wave observations (see 
e.g. the discussions by Testi et al. 2001; Wilner et al. 2005). 
To properly estimate the contribution of the gas emission to the 
long wavelength SED of CQ Tau, in our analysis we add to the 
dusty disk model a power law with variable spectral index ag 
as above, which we try to constrain using our new observations 
at long wavelengths with the assumption that all the emission 
we measure at 3.6 cm is due to gas emission. What we find 
with our analysis (which is presented in the next Section and 
in figures [2] and [3]l is only an upper limit to the possible gas 
contribution at 7 mm and 1.3 cm. 



4.3. Model fitting to the SED 

We used the dusty disk model and the parametrization for the 
gas emission discussed above to fit the submillimeter to cen- 
timeter wave SED of CQ Tau. In Table [2] we show the flux 
used at the various wavelengths from our own data and from 
the literature. The aim of SED fitting at mm-cm wavelengths 
is to constrain the maximum size of dust grains for the dust 
model we adopted (assuming a uniform dust population in the 
disk mid plane). To avoid degeneracies with the disk geometry, 
we assumed a disk outer radius of 63 AU, a disk inclination 
of inclination of 30°, and a surface density distribution index 
p = 0, as derived from the fit of the spatially resolved 1.3 mm 
dust emission as described in Section The inner radius is set 
to 0.2 AU, which corresponds to the dust evaporation radius as 



fitted minimizing the x as a function of the free parameters 
(a ma x, £ 40 , an d a s ). 

The left panel of Figure [2] shows the x 2 contours for aff 
and /3 assuming our reference dust. The data are consistent at 
the 68% confidence level with aff larger than -0.1 and /3 be- 
tween 0.49 and 0.68. For aff > 1 the free-free contribution at 
wavelengths up to 1.3 cm is negligible (less than ~IQ% at 1.3 
cm) and f3 becomes independent of ag. A gas emission index 
aff higher than -0.1 implies that the free-free contamination to 
the observed flux is at most 5% at 7 mm and 31% at 1.3 cm. 
Such a low value of the free-free contribution at 1.3 cm would 
be consistent with our non-detection of an unresolved emission 
due to the gas at at this wavelength (as shown in Section |2j), 
but, given the low signal to noise ratio of our 1.3 cm maps, fu- 
ture observations with higher sensitivity are needed to confirm 
this point. Our constraints on (3 are consistent and put on firmer 
grounds the previous estimates, which were 0.6±0.1 from |TestT| 



calculated in Isella et al. ( 200611. The disk atmosphere is consid- 
ered to be populated only with small dust grains with the same 
chemical composition as the midplane dust. The SED is then 



eTaLl ( [2003] l and 0.70±0.04 from Chapill on et al.| ( [2008| note 
that these latter authors compute the value of f3 between 1.3 
and 3 mm). Both of these previous studies could not constrain 
the cm emission and the free-free contamination at mm wave- 
lengths. 

As discussed by many authors (see e.g. the discussion in 
Natta et aT]2 007), while the constrain on the value of /3 is very 
strong from the millimeter SED fits, relating this to the exact 
properties of the dust grains in the disk is more difficult, due to 
the uncertainties in the dust composition (and hence opacity). 
Fig.[3]top left panel shows the constraint on the maximum grain 
sizes using our reference dust model, in this case we find that 
we can constrain the dust to be grown to ~1 cm pebbles. To 
explore how this conclusion depends on the dust properties we 
have also run a set of models varying the ice content, the level 
of porosity and the dust size distributions. In Fig. [3] we show 
the x 1 surfaces obtained exploring the constraints on a max and 
aff for a representative set of dust parameters. A summary of 
the results of all the fits is reported in Table|3] As expected, the 
grain sizes that we derive are larger in the case of more porous 
grains (no-ice and 50% vacuum cases) and are smaller for more 
compacted grains (10% vacuum case). Derived grain sizes are 
also smaller if we assume that the grain size distribution is more 
heavily weighted towards large grains (q=2.5 case); while for 
a size distribution with q= 3.5 very large grains, if present, do 
not contribute significantly to the dust opacity and we can only 
put a lower limit to the maximum grain size. The result of this 
analysis is that, even considering variations to the properties of 
the dust, our data is only consistent with maximum grain sizes 
larger than about 1 mm, implying that the grain growth process 
in this system is well advanced. In the remainder of this paper 
we will only consider the reference dust. 

5. Model fitting to the visibility data 

In the previous section we have shown that the grains on the 
CQ Tau disk midplane have reached maximum sizes of at least 
a few mm, confirming earlier suggestions. This result was ob- 
tained by assuming values of the disk outer radius, inclination 
and surface density index of 63 AU, 30° and 0, that are obtained 
as discussed below. 
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Table 3. Constraints on model parameters from the SED fitting. The range of values reported for each parameter is within the 
68% confidence level. 



Parameter 


ref. dust 


no ice 


vacuum 10% 


vacuum 50% 


q=2.5 


q=3.5 


a max (cm) 


0.8-2.4 


1.7-4.7 


0.1-0.3 


6.5-16 


0.2-0.7 


> 1 


as 


> -0.1 


> 0.0 


> -1.3 


> 0.1 


-1.1-1.2 


> -0.2 


X 40 (g cm" 2 ) 


2.6-5.2 


3.6-6.8 


0.55-1.05 


12.8-23.2 


1.3-2.7 


> 4.7 


f-fta. (%) 


<5 


<4 


£33 


<3 


1-24 


<5 


f-fl.3„„ (%) 


<31 


< 28 


< 100 


< 26 


9-84 


<35 



1 00.0 



10.0 



d8:o 



10.0 



10.0 



0. 1 



ref. dust 




vacuum; 



0% 



q = 2.5 




vacuum = 50% 




q = 3.5 




Fig. 3. Confidence limits on a max and as from the SED fitting 
using the different ensembles of dust; from top left to bottom 
right: reference dust with and without ice, volume percentage 
of vacuum 10% and 50%, grain size distribution index q-2.5 
and q-2>.5. Contours represent the 68% and 95% confidence 
levels for the fit. 
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the observed samples. Following the method described in Isella 
( 2007} , at each wavelength we computed the^- 2 function 



Following the procedure used in Testi et al. (2003 ), we pro- 
duced model images corresponding to different sets of param- 
eters: disk inclination and position angle (which we define in 
this paper as the angle east of north of the line defining the in- 
tersection between the disk plane and the plane of the sky), 
disk radius, surface density power law index p and normal- 
ization S40. For the dust properties we used the reference dust 
model discussed in the previous section for the disk mid plane 
and surface, assuming a constant composition along the entire 
disk. On the mid plane we adopted a maximum grain size of 
&max =1 cm ( as derived above). The model images of the disks 
at each wavelength were Fourier transformed and sampled at 
the appropriate positions on the (u,v) plane(s) corresponding to 



as: 

X 1 = Xj(\Rejiobs)-Rej{mod)f 

+ [lrrij(obs) — Imj(mod)] )/v 2 j, 

where Re and Im are the real and imaginary part of the observed 
(obs) and theoretical (mod) correlated flux, and <x is the noise 
on each measure. 

This procedure was repeated on a wide grid of model pa- 
rameters to construct a^ 2 hypercube. In order to find the best 
fitting model, the cube is then searched for the minimum of the 
X 1 as a function of all parameters. The procedure is repeated 
independently for all four datasets from 0.87 mm to 7 mm and 
the best fitting model is derived independently at each wave- 
length. This procedure has the advantage of allowing to search 
for systematic variations with wavelength of the model param- 
eters. The drawback is that at each wavelength we are fitting 
a large set of parameters and some of these are only weakly 
constrained by some of our datasets. The 1.3 cm dataset has a 
signal to noise ratio too low to perform a proper visibility fit- 
ting, while at 3.6 cm the emission is not due to the dust thermal 
emission in the disk, we thus excluded the two longest wave- 
length datasets from the visibility fitting procedure. 

In Fig. [4] we show the projections of the^ 2 hypercube on a 
few planes for each wavelength. The comparison between the 
observed images, the images reconstructed from the visibilities 
of the best fitting models, and the residual image obtained from 
the different between the observed and the best fitting visibili- 
ties are shown in Fig. [5] The comparison between the observed 
and best fitting real part of the visibilities is shown directly in 
Fig. [6] to produce the figure we deprojected the baseline lengths 
using the position angle and inclination of the disk as derived 
from the best fit, then we produced annular average visibilities 
for the data (points with error bars) and the model. In Table |4] 
we list the best fit values and the confidence range for the model 
parameters at each wavelength. 

The 1.3 mm dataset is the one with the best combination of 
resolution and signal to noise and is the one for which we derive 
a more constrained fit for the disk model parameters. The best 
fit values for the inclination and the position angle (measured 
east of north) are 34° and 42° respectively. The inclination is 
in good agreement with the value derived from observations of 
the CO emission (Chap illon et al.|2008 note that these authors 
report the position angle of the disk axis). The slope of the sur- 
face density distribution p is constrained to be between -0.2 and 
0, implying an almost constant radial dust distribution. Values 
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Fig. 4.x 2 hypercube projections on the (p,Logio(X4o)), (p,R ou tX 
and (Logio(S4o),R ou t) planes, to produce each plot we have cho- 
sen the value of the other parameters corresponding to the min- 
imum of x 1 at each position on the planes shown. Contours rep- 
resent the 68% and 95% confidence levels for the fit. In each 
row we show the results for the fit at a given wavelength, from 
top to bottom ar 0.87 mm, 1.3 mm, 2.7 mm and 7 mm (as la- 
belled to the right). 

Table 4. Results of the model fits to the visibilities. For each 
parameter we report the best fit value and the minimum and 
maximum values of the 68% confidence level intervals (lcr). 



Parameter 



£ 4 0AU (g cm ) 
P 

Rout (AU) 



0.87 mm 



1.3 mm 



2.7 mm 7 mm 



7.1 + 



0.04!™ 



5A + _f 



2.5 + 



2.5 + 



-0 075+ 0025 5 +0 - 3 

u.u/j 0125 v.j^ 15 



63 + 



70+i° 

— o 



5 +3) 

J - -3.0 

0.5 +0 - 5 



52 +33 



of p close to are also obtained from the observations at 0.87 
mm, while p ~ 0.5 is suggested by the observations at longer 
wavelengths. As discussed in Isella et al. ( 2010| l, a systematic 
variation of the fitted surface density distribution from 0.87 and 
7 mm would suggest a radial variation of the slope of the dust 
opacity /3. Values of p increasing with the wavelength would 
imply that /3 increases with the distance from the central star 
(from Eq. 9 in Isella et aI7||2.Q10] >, or that the maximum grain 
size decreases with the orbital distance. In particular, a varia- 



< 




residuals 




residu&ls"""' 
■T l 1 1 I i l l I fT) I i I h 



residuals. 



Aa (") 

Fig. 5. Comparison between the observed images at 0.87 mm, 
1.3 mm, 2.7 mm and 7 mm (left panels, from top to bottom) and 
the images corresponding to the best fit model at each wave- 
length obtained by sampling the model visibilities at the same 
locations and with the same noise as for the observed data (mid- 
dle panels). The images derived by the difference from the ob- 
served and model visibilities are shown on the right panels. At 
each wavelength the contours are the same in all three plots 
and correspond to -2cr and from 2<x to the peak value spaced 
by 2<x, where <x is the rms measured in the observed map (see 
also Fig.[T|. 

tion of p of about 0.5 between 0.87 and 7 mm, corresponds to 
an increase of (3 from to about 1.4, between 0.2 and 63 AU. 
In terms of grain size distribution, this would imply that the in- 
ner disk is mainly populated by dust grains of few centimeters 
in radius, while the outer disk is characterized by a sub-micron 
grain size distribution similar to what observed in the inter- 
stellar medium. Our results are in qualitative agreement with 



recent model predictions by |Birnstiel et al. (2010ai who pre- 
dict values of (3 in circumstellar disks to grow as a function of 
radius due to the different growth level of grains as a function 



of radius as predicted by the dust evolution models of Brauer 



|et al. | ( 12008) and |Birnstiel et al.| ( |20T0bl >. Nevertheless, our ini- 
tial results should be taken with caution: with the exception of 
the 1.3 mm observations, the angular resolution and the signal 
to noise level is too low for us to constrain reliably the value 
of p at the various wavelengths, and more sensitive high angu- 
lar resolution observations together with a proper fitting with 
a disk model which takes into account the dust variations as a 
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Fig. 6. Deprojected radial averages of the real part of the visi- 
bility as a function of baseline length. The full lines show the 
predicted visibilities for the best fit models at each wavelength. 



function of radius are necessary to confirm our tentative con- 
clusions. 

From our observations, we measure an outer disk radius of 
about 60 AU, which is much smaller than the value of 140±14 
AU derived from CO and dust continuum emission observa- 
tions at 2" resolution (Chapillon et al. 2008 note that we scaled 
the values of their Table 3 for a stellar distance of 100 pc). 
When using a simple power-law description for the surface 
density, it is well known that disk radii mesured from the dust 
millimeter continuum emission are often smaller than those de- 



rived from molecular line observations (e.g. Hughes et al. 2008 



Isell a et al.|2007| l. Consequently, the inferred radius of 60 AU 
has to be considered as a lower limit. 



6. Summary and conclusions 

In this paper we have presented new high angular resolution ob- 
servations from the submillimeter through the centimeter wave- 
length range of the CQ Tau system. At millimeter and submil- 
limeter wavelegths we detect and resolve the thermal dust emis- 
sion from the circumstellar disk. At 3.6 cm we detect compact 
emission from the inner regions of the system that we attribute 
to gas emission. We fit the observations as the emission of a 
dusty disk with an additional contribution of unresolved gas 
emission. The combination of high angular resolution observa- 
tions at several wavelengths and the accurate estimate of the 
possible contamination due to gas allows us to derive more ac- 
curately the properties of the thermal dust emission from the 
circumstellar disk itself. 

We find that the average value across the disk of the dust 
opacity coefficient computed between 1.3 and 7 mm, is 
~ 0.6 ±0.1. This confirms and puts on firmer ground previous 
derivations of [3 in this system. These values of (3 imply that sig- 
nificant grain growth has occurred in this disk, with the largest 
grains grown to ~1 cm. The exact value of the maximum grain 



size depends on the details of the dust model adopted and we 
find it likely to be in the range from 1 mm to 10 cm. 

Our resolved images at several (sub)millimeter wave- 
lengths show a weak indication of the signature of grain growth 
variation as a function of radius in the CQ Tau disk. We thus 
find marginal evidence to support the expectation that the in- 
ner regions of the disk contain grains grown to larger sizes than 
the outer regions. This result needs to be confirmed by higher 
signal to noise data at the extremes of the wavelength range 
probed by our observations. In particular we do expect that in 
the near future EVLA observations at 7 mm and 1.3 cm will 
allow to confirm (or disprove) our suggestion. 
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